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ABSTRACT 

We construct merger trees of dark matter haloes and quantify their merger rates and 
mass growth rates using the joint dataset from the Millenniiim and Millcnni\im-II 
simulations. The finer resolution of the Millennium-II Simulation has allowed us to 
extend our earlier analysis of halo merger statistics to an unprecedentedly wide range of 
descendant halo mass (10^° < Mo < IQI^Mq), progenitor mass ratio (IQ-^ < ^ < 1), 
and redshift (0 ^ z < 15). We update our earlier fitting form for the mean merger 
rate per halo as a function of Mq, ^, and z. The overall behavior of this quantity 
is unchanged: the rate per unit redshift is nearly independent of z out to z ~ 15; 
the dependence on halo mass is weak (oc Mq-^^); and it is nearly a power law in the 
progenitor mass ratio (oc ^~^). We also present a simple and accurate fitting formula 
for the mean mass growth rate of haloes as a function of mass and redshift. This 
mean rate is AQMq yr^^ for IQ^'^Mq haloes at z = 0, and it increases with mass as 
oc M^-^ and with redshift as (1 + z)^-^ (for z > 1). When the fit for the mean mass 
growth rate is integrated over a halo's history, we find excellent match to the mean 
mass assembly histories of the simulated haloes. By combining merger rates and mass 
assembly histories, we present results for the number of mergers over a halo's history 
and the statistics of the redshift of the last major merger. 



1 INTRODUCTION 

Mergers of dark matter haloes are intimately connected to 
a wide array of phenomena in the now-standard ACDM 
cosmology. In addition to being the dominant channel for 

mass growth of haloes themselves, mergers are also respon- 
sible for the growth of stellar mass in galaxies, both directly 
via galaxy-galaxy mergers, and indirectly via the accretion 
of potentially star-forming gas. Furthermore, mergers help 
shape many important observational properties of galaxies, 
e.g., star formation rates, color and morphology transforma- 
tions, dynamical states of stellar disks, and galaxy mass and 
luminosity functions. Mergers are also responsible for the ex- 
istence of satellite galaxies such as dwarf spheroidals in the 
Milky Way and non-cD galaxies in galaxy clusters. Quanti- 
fying the rate of halo-halo mergers, and its possible depen- 
dence on factors such as halo mass, mass ratio, and time, is 
therefore of great interest for a theoretical understanding of 
galaxy formation and its connections to observations. 

In a series of papers, we have examined various as- 
pects of the growth of dark matter haloes. In Fakhouri & 
Ma (2008), we computed the merger rates of dark mat- 
ter haloes from the Millennium Simulation (Springel et al. 
2005) and presented a simple algebraic fitting form for our 
results. The resolution and size of this simulation allowed 
us to determine the merger rate over the parameter range of 
10^^ < Mo < lO^^M© for the mass of the descendant haloes, 
10~^ J; C ^ 1 for the mass ratio of the progenitor haloes. 



and ^ z < 6 for the redshift. The detailed environmental 
dependence of the merger rates and halo mass growths was 
analyzed in two subsequent papers (Fakhouri & Ma 2009, 
2010). In McBridc ct al. (2009), wc studied the statistics 
of the halo mass assembly histories and mass growth rates 
in the Millennium Simulation. Halo mergers have also been 
studied in a handful of papers by others (e.g., Governato 
et al. 1999; Gottlober et al. 2001; Berrier et al. 2006; Mailer 
et al. 2006; Guo & White 2008; Genol et al. 2009; Stewart 
et al. 2009). The pre-2008 studies were all limited to small 
simulations that mainly investigated major mergers in a nar- 
row mass range at low redshift (typically z < 1). Some such 
studies have emphasized potential challenges for hierarchical 
structure formation; for instance, Stewart et al. (2008) have 
noted that the frequency of major mergers among Milky- 
Way sized haloes poses a problem for thin-disk survivabil- 
ity. Much work has also been done in quantifying halo mass 
accretion and assembly histories using A'^-body simulations 
that are smaller than the Millennium runs (e.g., Laccy & 
Cole 1994; Tormen et al. 1997; Tormen 1998; Wechsler et al. 
2002; van den Bosch 2002; Li et al. 2007; Zhao et al. 2009, 
except Cole et al. 2008). 

In this paper, we extend the results presented in 
Fakhouri & Ma (2008) and McBride et al. (2009) by in- 
corporating the Millennium-II Simulation (Boylan-Kolchin 
et al. 2009). This simulation has the same number of par- 
ticles as the Millennium Simulation but has 125 times bet- 
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ter mass resolution. This new database provides 7.5 x 10 
dark matter haloes (each containing more than 1000 simula- 
tion particles) between redshift and 15 and their subhalo 
merger trees for our analysis. Adding to the 11.3 x 10^ haloes 
(between z = and 6) available from the Millennium Sim- 
ulation, this combined dataset allows us to determine the 
dark matter halo merger rates and mass growth rates from 
« = to up to z = 15, for over five orders of magnitude 
in the descendant halo mass (10^° < Mo < lO^^Mo) and 
progenitor mass ratio (10~* ^ ^ ^ !)• 

This paper is organized as follows. Section 2 describes 
the dark matter haloes in the Millennium and Millennium- 
II simulations, and how we construct the merger trees and 
quantify the merger statistics and mass accretion histories 
of the haloes. In Section 3, we present results for three types 
of statistics: merger rates at 2 = up to ~ 15 for halo mass 
~ 10^° to IO^'^Mq (§ 3.1); the rate at which the haloes are 
accreting dark matter across the virial radii, and the mass 
growth history of haloes (§ 3.2); and the cumulative merger 
statistics over a halo's past history, e.g., the mean cumula- 
tive number of mergers of a given mass ratio experienced as 
a function of z and halo mass, and the distribution of the 
redshift at which the last major merger occurred for haloes 
at various mass and redshift (§ 3.3). The Appendix contains 
a detailed comparison of the three types of algorithms that 
we have tested for handling the fragmentation events in a 
merger tree of FOF haloes (Fakhouri & Ma 2008, 2010). 
They are named "snip," "stitch," and "split," depending on 
whether the fragmented subhalo was ignored, stitched back 
to the original FOF halo in subsequent outputs, or split off 
from the FOF at earlier times. A quantitative assessment of 
the systematic differences in the merger rates derived from 
each algorithm is provided in the Appendix. 

The cosmology used throughout this paper is identical 
to that used in the the Millennium simulations: a ACDM 
model with D,m = 0.25, fit = 0.045, D.a = 0.75, h = 0.73, 
an initial power-law index n = 1, and ag — 0.9. Masses and 
lengths are quoted in units of Mq and Mpc without the 
Hubble parameter h. 



2 CONSTRUCTION OF HALO MERGER 
TREES 

2.1 The Two Millennium Simulations 

The Millennium and Millennium-II simulations are large N- 

body simulations of cosmological structure formation using 
the concordance ACDM cosmological parameters listed at 
the end of Section 1. The simulations are described in detail 
in Springel et al. (2005) (Millennium) and Boylan-Kolchin 
et al. (2009) (Millennium-II); here we summarize some basic 
features of the simulations and of the default post-processing 
procedures that result in subhalo merger trees. 

Both simulations follow the evolution of 2160'"* ~ 10^" 
particles from redshift 127 to redshift using versions of the 
GADGET tree-PM code (Springel et al. 2001; Springel 2005). 
The sirrmlations differ in spatial scale and mass resolution: 
the Millennium Simulation uses a box size oi L = 685 Mpc 
and a Plummer-equivalent force softening that is a fac- 
tor of 10^ smaller, e = 6.85 kpc, with a particle mass of 
rup = 1.18 X 10^ Mq. The Millennium-II Simulation uses 



L = 137 Mpc and e = 1.37 kpc, both of which are a fac- 
tor of 5 smaller than the values from the Millennium Sim- 
ulation; the particle mass is therefore 125 times smaller, 
TJip = 9.43 X lO*' Mq . The two simulations have 60 outputs 
at identical redshifts between « « 20 and 2 = 0, spaced ap- 
proximately equally in log 2, as well as additional snapshots 
(4 for the Millennium, 8 for the Millennium-II) at higher 
redshifts. 

Subhalo merger trees are constructed in an identical 

fashion for the Millennium and Millennium-II simulations. 
Dark matter haloes are first identified at each snapshot using 
a Friends-of-Fricnds group-finder (FOF; Davis et al. 1985) 
with a linking length of 0.2 times the mean interparticle sep- 
aration. All FOF groups with at least 20 particles are stored. 
The SUBFIND algorithm (Springel et al. 2001) is then applied 
to each FOF group to identify halo substructure. SUBFIND 
identifies local density maxima and performs an unbinding 
procedure to determine which particles in the FOF group 
are bound to each density peak. Substructures with at least 
20 particles after unbinding arc stored, resulting in a list of 
subhaloes (SHs) associated with each FOF group in the sim- 
ulation. Note that some FOF groups do not contain 20 self- 
bound particles and therefore not every FOF group contains 
a subhalo, while some FOF groups can contain many self- 
bound density peaks and therefore have many subhaloes. 

These subhaloes are then linked across simulation snap- 
shots to produce subhalo merger trees. This linking is done 
by establishing a unique descendant for each subhalo in the 
following manner. First, all particles in the subhalo are rank- 
ordered by binding energy. A list of candidate descendants - 
all subhaloes at the subsequent snapshot containing at least 
one particle from the subhalo in question - is built and a 
figure of merit is computed for each descendant. The can- 
didate descendant with the highest score is assigned as the 
actual descendant. The figure of merit for candidate descen- 
dants is simply a weighted sum of the rank-ordering of the 
subhalo's particles; this procedure ensures that the tightly 
bound center of a subhalo is weighted more heavily than the 
less-bound outer regions even if the center is subdominant 
in terms of mass. 

In addition to searching for a descendant at the subse- 
quent output, a search is also performed two snapshots later. 
This additional step accounts for subhaloes that are tem- 
porarily unresolved when passing near the center of a more 
massive system but re-appear later. On occasion, no descen- 
dant can be identified at either of the two subsequent snap- 
shots, in which case, the subhalo is not assigned a descendant 
at all but rather is considered destroyed. With subhaloes and 
their unique descendants identified, subhalo merger trees are 
built by linking subhaloes and their descendants: all sub- 
haloes with a common descendant at 2 = are linked to 
all subhaloes sharing these subhaloes as descendants, and 
so on. A given subhalo merger tree thus contains all sub- 
haloes that can be linked via their descendants to one spe- 
cific subhalo at 2 = 0. The trees link 760 million subhaloes 
for the Millennium Simulation and 590 million subhaloes 
for the Millennium-II Simulation. For the central subhalo 
of a 2 = galaxy-mass halo (M « 10^^ M©), its subhalo 
merger tree typically consists of 90 subhaloes in the Millen- 
nium Simulation and 2800 subhaloes in the Millennium-II 
Simulation. 
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Table 1. The number of merger events in the two Millennium simulations at four representative redshifts {z si 0, 0.5, 1, and 2). At each 

2, wc list the number of descendant FOF haloes that have a single progenitor halo {Np = 1, i.e., no mergers) and multiple progenitors 
(Np ^ 2), for five descendant mass bins (left to right). The descendant mass here refers to the halo mass at the redshift listed rather 
than at the present day. Only haloes containing more particles than our minimum cutoff (1000 for descendants; 40 for progenitors) are 
counted. The higher-resolution Millennium-II Simulation dominates the contribution to the merger statistics of Mq < 10^^ Mq haloes, 
while the larger-volume Millennium Simulation dominates the contribution to cluster-mass haloes. 



2.2 Halo Fragmentation 

In this paper, as in our previous work (Fakhouri & Ma 2008; 
McBride et al. 2009) , our focus is on the merger and assem- 
bly histories of FOF haloes. To do this we must first con- 
struct merger trees of the FOF haloes from the underlying 
subhalo trees described in Section 2.1. Such construction is 
nontrivial due to halo fragmentations: subhalocs of a pro- 
genitor FOF halo may have descendants that reside in more 
than one FOF halo. Sometimes this is duo to a physical un- 
binding event in which a subhalo formerly bound to an FOF 
is ejected out of the FOF system. Sometimes the fragmenta- 
tion is spurious a subhalo may oscillate in and out of the 
FOF group before finally settling in. Sometimes the FOF al- 
gorithm incorrectly groups subhaloes that are unbound but 
only happen to pass by one another and should not be as- 
sociated as a single FOF group. 

Wc presented detailed comparisons in Fakhouri & Ma 
(2008, 2010) of three types of algorithms - snip, stitch, and 
split - for handling these fragmentation events. In the Ap- 
pendix we summarize these algorithms and quantify the sys- 
tematic differences in the merger rates derived from each 
algorithm. 

For the main results presented in § 3 below, we use the 
split-3 algorithm, in which the subhalo fragments that pop 
out of an FOF halo are either snipped or split depending on 
a simple criterion. The fragmented subhalo is snipped if it 
is observed to remain in the FOF halo for all 3 snapshots 
immediately preceding the fragmentation event; in this case 
the ancestral link between fragment and FOF is severed. If 
the fragmented subhalo is not in the FOF halo for all 3 pre- 
ceding snapshots, it is interpreted as distinct and is split off 
from the FOF. The split-3 algorithm generally gives very 
similar results to the stitch-3 algorithm used in Fakhouri & 
Ma (2008), e.g., the two methods produce merger rates that 
agree to within 10% for the redshifts and mass ranges that 
we have statistics for. The only exception is in the major 
merger regime for low-mass haloes at low redshift {z < 1), 
where split-3 is lower than stitch-3 by up to 30% (sec Fig. 10 
in Appendix). Overall, split-3 appears slightly more robust 
at handling spurious FOF linking events in this regime (also 
see Fakhouri & Ma 2010; Genel et al. 2009). As discussed in 
the Appendix, however, the exact definition of what consti- 



tutes a merger may be situation-dependent, meaning that 
no single method is perfect in all cases. 



2.3 Extracting Merger Rates and Mass Accretion 
Histories 

From the merger trees of FOF haloes obtained by applying 
a given fragmentation algorithm, we extract a merger cata- 
log. Each catalog provides us with a list of descendant FOF 
haloes at redshift Zd ^ with mass Mq , and for each descen- 
dant halo, its set of Np FOF progenitors at Zp — Zd + Az, 
where Np can range from 1 (i.e. a single progenitor) to a 
large number, depending on the halo mass and the value of 
Az. We label the rank-ordered progenitor mass with M,, 
i G (1,2, . . . ,iVp), and Mi > A/2 > ■ ■ ■ Mn^- To ensure 
that only numerically resolved haloes are included in our 
study, we impose a minimum of 1000 particles for the de- 
scendant haloes and 40 particles for the progenitor haloes. 
For the Millennium Simulation, this criterion corresponds to 
a minirrmm halo mass of 1.2 x 10^2 for the descendant 
and 4.7 X lO^" Mq for the progenitor. For Millennium-II, the 
minimum masses are 9.4 x 10® M0 and 3.8 x 10* M0 for the 
descendant and progenitor haloes. We emphasize that the 
mass of a descendant halo refers to its mass at a given red- 
shift Zd and not its ultimate mass at 2 = (unless Zd ~ 0). 

We compute the merger rates at redshift z as a func- 
tion of descendant mass Mo and progenitor mass ratio 
C = M^/Ml (for i > 1). We define B{Mo,^,z) to be the 
number of mergers per Mpc'^, dMo, d^, andAz with mass 
Mo ± dMo/2 and mass ratio ^ ± As discussed in 

Fakhouri & Ma (2008), we find the mean merger rate per 
halo, B{Mo,^,z)/n{Mo,z) = dNm/d^/dz, where n{Mo,z) 
is the rmmbcr density of haloes, to have a particularly sim- 
ple dependence on the merger parameters. This rate, when 
expressed in per redshift units, is a dimensionless quantity 
that gives the mean number of mergers per halo per unit z 
per unit ^. To avoid artificial boundary effects at « = 0, we 
use the two outputs at z = 0.12 and 0.06 to compute the 
a ~ merger rate. 

To compute the mass accretion histories and accretion 
rates of haloes, we start with a given descendant FOF halo 
at some redshift and identify the mass of its most massive 
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progenitor at an earlier snapshot. This process is iterated 
backwards in time to construct the main branch of the de- 
scendant's merger tree. The mass trajectory along the main 
branch of a descendant gives us its mass accretion history 
M{z) (see. e.g., Lacey & Cole 1993), from which we can 
compute the mass accretion rate M as a function of z. Note 
that the progenitor halo on the main branch of a descen- 
dant halo at a given snapshot need not be the most massive 
progenitor of that descendant at that snapshot. 

2.4 Definitions of Halo Mass 

In our prior analysis of the Millennium halo merger rate 
(Fakhouri & Ma 2008), we assigned the halo mass using 
the standard FOF mass Mfof- This mass is simply pro- 
portional to the number of particles assigned to each FOF 
halo by the FOF group finder. An alternative definition 
that we will use throughout this paper is the sum of the 
masses of an FOF's subhaloes, Msh- This definition has 
been shown recently by Genel et al. (2009) to be more ro- 
bust than Mfof since the SUBFIND algorithm assigns only 
gravitationally bound particles to each subhalo. 

Overall, we find the halo mass functions computed using 
these two mass definitions to differ at the 5% level at all 
halo masses. This difference can be caused by a shght excess 
of mass in FOF haloes due to unbound or spuriously linked 
particles, as well as by a slight deficit in Msh when SUBFIND 
does not account for all the mass physically associated with 
a subhalo. When restricted to the subset of haloes that are 
undergoing very minor mergers, however, Genel et al. (2009) 
noted that the FOF mass of the smaller progenitor increases 
as it approaches the more massive progenitor. For minor 
mergers involving mass ratios as low as ^ ~ 0.001, the ratio 
Mfof /Msh for the smaller progenitor can rise from 1.03 
to 1.5 prior to mergers. We will therefore use Msh for halo 
masses in this study. We note that this discrepancy occurs 
only for the small subset of low-mass haloes that are in the 
process of merging onto a much larger halo; its effect on the 
total halo mass function is therefore limited to ~ 5%. 



3 RESULTS 

To provide a sense for the halo statistics and merger events 
available from the two Millennium simulations, wo list in 
Table 1 the number of descendant haloes (above 1000 parti- 
cles) and their progenitors (above 40 particles) at four rep- 
resentative redshifts for five broad mass ranges from 10^° to 
> 1O^*M0. The results presented below are based on these 
events and those at other redshifts. 

3.1 Merger Rates 

3.1.1 Present-Day Merger Rates 

The left panel of Fig. 1 presents B{Mo,$,, z = 0), the z = 
mean number of mergers per unit volume, descedant mass 
Mo, mass ratio^, and redshift as a function of progenitor 
mass ratio £, from the two Millennium simulations (solid for 
Millennium II; dashed for Millennium). The colored curves 
correspond to different mass bins ranging from 10^° M© 
(blue) to 10^® Mo (red). The rates are determined from the 



z = 0.06 and 0.12 merger tree catalogue since, as described 
in Section 2.3, we would like to avoid the z = 0.0 snapshot 
due to the boundary effects that interfere with the post- 
processing algorithms used to handle the halo fragmentation 
events. The split-3 algorithm is used here; other algorithms 
yield qualitatively similar agreement between the two simu- 
lations (see Appendix for details). 

The right panel of Fig. 1 shows the mean merger rate 
per halo, B/n = dNm/ d^/ dz, where each of the curves in the 
left panel has been divided by the number density of haloes 
in that mass bin. The collapse of the curves to nearly a sin- 
gle curve shows that the per halo merger rate dNm/d(,/dz 
is nearly independent of halo mass. This collapse is similar 
to that seen in Fig. 6 of Fakhouri & Ma (2008) for the Mil- 
lennium Simulation. A comparison of the two figures helps 
illustrate the large dynamic range achieved when the two 
Millennium simulations are combined: the halo mass range 
has been increased by two orders of magnitude in Fig. 1, 
and for each mass bin, the progenitor mass ratio is extended 
downward by also a factor of ~ 100, reaching ^ ~ 10~^ for 
Mo = W^^Mq. 

The overlap in the merger parameter space between the 
two simulations is seen to be fairly small in Fig. 1. The 
two simulations are therefore quite complementary: Millen- 
nium II allows us to probe descendant and progenitor masses 
that are a factor of 125 smaller than Millennium, whereas 
the larger box of the Millennium Simulation provides robust 
statistics for the rare events that are poorly sampled in Mil- 
lennium II, e.g., major mergers of massive haloes (i.e. the 
lower right corner of left panel of Fig. 1). Over the small re- 
gion of overlap. Fig. 1 shows good agreements between the 
merger rates determined from the two simulations: both the 
power-law dependence on ^ and the weak dependence on 
Mo carry over from Millennium to Millennium II. Boylan- 
Kolchin et al. (2009) show that many other quantities, such 
as halo mass functions, formation times, and subhalo abun- 
dances, have a much wider range of overlap and that the two 
simulations are in excellent agreement for these quantities 
as well. 

The weak dependence of the merger rate on Mo is shown 
explicitly in Fig. 2. Each curve here shows the mean rate 
per halo, dNm/d^/dz, integrated over different ranges of 
^ ^ ^min, where ^mm = 0.3,0.10.01, and 10~^ (from bot- 
tom up). Major mergers with mass ratio within 1:3 (bottom 
curve) are clearly much more rare than minor mergers (top 
curves), but all the curves have very similar power-law de- 
pendence on Mo. Over about 4.5 orders of magnitude in Mq, 
the rate increases by only a factor of ~ 3, suggesting that 
the merger rate scales roughly as ~ Mq'^. A more accurate 
fit is provided in Sec 3.1.3 below. 

3.1.2 z> Merger Rates 

The Millennium Simulation provided sufficient halo statis- 
tics for us to determine the halo merger rates up to « ~ 6 
in our previous study. The higher mass resolution of the 
Millennium-II Simulation now allows us to probe redshifts 

up to ~ 15. The combined results from the two simulations 
arc shown in Fig. 3, which plots the mean merger rate per 
unit redshift (left panel), dNm/dz, and per unit time (right 
panel), dNm/dt, as a function of redshift. These merger rates 
have been integrated over different ranges of ^ ^min, rang- 
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Figure 1. Left panel: The mean merger rate of 2: = FOF haloes, B{Mo,^), as a function of the mass ratio of the progenitors (§) 
and the descendant halo mass (Mq) over 6 orders of magnitude: 10^" to IO^^^q from right (blue) to left (red). The Millcnnium-II 
results are shown in solid, while the results from the Millennium are in dashed curves. Right panel: The mean merger rate per halo, 
B{Mo, ^) /n{Mo) = dNm/d^/dz. Normalizing B(Afo,C) by the halo number density n(Mo) collapses the curves in the left panel to nearly 
a single curve, indicating that dNm/d^/dz is nearly independent of Mg and has a simple universal form. 
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Figure 2. The z ^ mean merger rate per halo (per unit z), 
dNm/dz, as a function of descendant mass, Mq, for four ranges 
of progenitor mass ratio ^. The Millennium-II results are shown in 
solid, while the original Millennium results are the dashed curves. 
The upper curves include increasingly more minor mergers. The 
mass dependence is weak over five orders of magnitude in mass 
and is well approximated by a power law oc M^-^33 



ing from major mergers with ^min = 0.3 (solid curves at 
bottom), to extreme minor mergers with ^min ~ (dot- 
ted curve at top). Within each line type, the colors indicate 
different descendant mass bins ranging from 10^'' (blue) to 
> 1O"M0 (red). Only the higher mass bins are plotted as 
Cmin is lowered. This is because minor mergers of low-mass 
haloes fall below the mass resolution limit. 

Fig. 3 indicates that the general trends reported in 



Fig. 8 of Fakhouri & Ma (2008) continue to hold in the 
Millennium-II Simulation. The dimensionless rate dNm/dz 
is remarkably independent of redshift up to z ~ 15, whereas 
the rate per Gyr, dNm/dt, rises with increasing z because 
a unit redshift corresponds to a shorter time interval at 
higher z. This redshift dependence is similar to that ob- 
tained by Guo & White (2008) for the dimensionless growth 
rates due to mergers of both haloes and galaxies (based on 
semi-analytic models) from the Millennium Simulation. 



3.1.3 Merger Rate Fitting Forms 

Since the merger statistics in Figs. 1-3 are very consistent be- 
tween the two Millennium simulations, we use an analytical 
form similar to equation (12) of Fakhouri & Ma (2008) to fit 
the dimensionless mean merger rate dNm/ d^/ dz (in units of 
mergers per halo per unit redshift per unit ^) from the com- 
bined Millennium dataset. An appealingly simple feature of 
this fitting form is that it is separable in the variables Mo, 
and z: 



dNrr 



{M,^,z) = A 



M 



t, exp 



(1) 

We find the best-fit parameters to be (a, /3, 7, rj) — 
(0.133,-1.995,0.263,0.0993) and {A,i) = (0.0104, 9.72 x 
10-3). The near z-independence in the left panel of Fig. 3 is 
more striking than in our 2008 study due to the larger cover- 
age in redshift here. In view of this lack of z-dependence, we 
choose to use the simpler factor of (1 + z)^ here rather than 
the growth rate of the density field used in Fakhouri & Ma 
(2008). In comparison to our 2008 study, the power-law slope 
of the mass dependence is slightly steeper here (a — 0.133 
vs. 0.089), whereas the power-law slope of the ^ dependence 
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Figure 3. Time evolution of the mean halo merger rates per halo in units of per redshift, dNm/dz (left panel), and in units of per Gyr, 
dNm/dt (right panel) from the two Millennium simulations. The descendant mass Mq and progenitor mass ratio ^ over five orders of 
magnitude are plotted. The weak dependence of the rates on Mq is shown by the different colors: 10^'' (blue), 10^^ (cyan), 10^^ (green), 
lO'^^ (orange), to > IO'^'^Mq (red). The line types denote different types of mergers, ranging from major mergers (solid) to extreme minor 
mergers (dotted) The rate dNm/dz on the left is remarkably constant out to ^ ~ 15; the rapid rise of dNm/dt with increasing z on the 
right is therefore largely due to the cosmological factor dt/dz, which spans a shorter time per unit z with increasing z. 



is slightly shallower here (/3 = -1.995 vs. -2.17). These dif- 
ferences are primarily due to the differing definitions of halo 
mass used in the two studies (FOF vs. sum of subhalos; see 
Sec. 2.4) and the refinements in our stitch-3 algorithm (see 
the Appendix). 

We note that the left panel of Fig. 3 does show mild 
variations in the redshift dependence among the different ^ 
bins: the rate increases slightly with increasing z for ma- 
jor mergers, while it declines somewhat for the very minor 
mergers (^min ~ 10"** to 10~^). Since this variation is so mi- 
nor and the minor merger regime is more prone to numerical 
resolution issues, we have opted for simplicity rather than a 
more complicated fitting form. 



3.2 Mass Growth Rates and Assembly Histories 

In the last section, we presented results for the instantaneous 
rates of halo mergers as a function of redshift, descendant 
mass, and progenitor mass ratio. Here, we examine a related 
set of statistics that quantify the mass growth of haloes. 
These two quantities are clearly related since mergers are 
a primary process for haloes to gain mass, but mergers are 
not the only process. As discussed at length in Fakhouri & 
Ma (2010), "diffuse" accretion of unresolved haloes or dark 
matter particles also makes an important contribution to 
halo growth. In the mass assembly history of a halo, mergers 
with other haloes typically result in more discrete but less 



frequent changes in the halo mass, while diffuse accretion 
leads to a more continuous change. 



3.2.1 Mass Accretion Rates 

To compute the total mass growth rate of a halo of a given 
mass Mo at time t, we follow the main branch of its merger 
tree (see § 2.3) and set M = (Mq - Mi) /At, where Mq is 
the descendant mass at time t and Mi is the mass of its 
most massive progenitor at time t — At. The mean value of 
M as a function of z for the complete set of resolved haloes 
in the two Millennium simulations is plotted in Fig. 4 (solid 
curves) . Nine ranges of Mo spanning five orders of magnitude 
(10^° Mq to IO^'^Mq from bottom up) are shown. Fig. 4 
can be compared directly to Fig. 5 of McBride et al. (2009) 
for the Millennium Simulation alone. The rising (AI) with 
increasing redshift in our earlier study is seen to continue to 
z ~ 14, and the nearly linear scaling of (M) with halo mass 
is extended down to ~ 10^°Mq. 

We find the mass accretion rates shown in Fig. 4 to be 
very well fit by the forms given by equations (8) and (9) 
of McBride et al. (2009). The coefficients quoted there only 
need minor adjustments after the Millennium-II results are 
added. We suggest the following updated fits to the mean 
and median mass growth rates of haloes of mass M at red- 
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Figure 5. Differential (top) and cumulative (bottom) distributions of the baryonic accretion rates, Mi,, for halo masses 10^" (left), 10^^ 
(middle), and > 10^*Mq (right). Within each panel, the accretion rates a,t z = 0,0.5, 1, and 2 are shown (except the right panel, where 
such massive haloes are present only at z = and 0.5). The distributions are seen to broaden significantly with increasing z. The vertical 
axis in each of the bottom panels labels the number of haloes per comoving Mpc^ that are accreting at a rate of Mi, or above. 
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Figure 4. Mean mass accretion rate of dark matter onto haloes 
as a function of redshift from the two Millennium simulations 
(solid curves). Halo masses ranging from IO^^A/q to > IO^^A^q 
are plotted. The dashed curves show the accurate fit provided 
by equation (2). The right-hand side of the vertical axis labels 
the mean accretion rate of baryons. Mi,, assuming a cosmological 
baryon-to-dark matter ratio of 1/6. 

shift z: 

= 46.1Meyr-(^^)" 

x{i + i.iiz)x/a4iTJpTnI 

x{l + 1.65z)^n,^{l + 2)3 + Oa . (2) 



At a given mass and redshift, the mean rate is overall higher 
than the median rate since the distribution of M has a long 
positive tail (see Fig. 5). The dashed curves in Fig. 4 illus- 
trate the remarkable accuracy of this formula in matching 
the simulation results over the broad ranges of halo mass 
and redshift shown. 

The right-hand-side label along the vertical axis of 
Fig. 4 shows the corresponding mean accretion rate of 
baryons, Mb, where we have assumed a cosmological baryon- 
to-dark matter ratio of ^It/^m = 1/6. These values are 
meant to provide a rough approximation for the mean rate 
at which baryons are being accreted near the virial radius 
of a dark matter halo. Most of these baryons are presum- 
ably in the form of warm or hot ionized hydrogen gas that 
is being channeled into the haloes along cosmic filaments, 
and various gas cooling and feedback processes will likely 
affect the baryon accretion rate. Many studies on galaxy 
formation are aimed at quantifying these physical processes 
under which these baryons are cooled to form neutral gas, 
molecular gas, and stars, and the feedback processes that 
heat up the baryons and lead to large-scale outflows. 

In Fig. 5 we show the differential (top) and cumula- 
tive (bottom) distributions of the baryonic accretion rate 
for three halo masses (left to right panels) and four red- 
shifts. A cosmic ratio of ilb/^m = 1/6 is again assumed to 
convert the dark matter rate into a baryonic rate. The distri- 
butions are strongly peaked at the mean values presented in 
Fig. 4 but exhibit long tails towards high positive values due 
to major merger events and towards negative values due to 
tidal stripping and halo fragmentation. Not only is the mean 
accretion rate higher at higher z, the distribution of Mb is 
also broader at higher z. For example, the comoving den- 
sity of Milky Way-mass haloes that are accreting baryons 
at a rate of at least lOOM© per year is approximately ten 
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Figure 6. Top panel: Mean mass assembly history M{z) of all 
2 = resolved dark matter haloes in the two Millennium simula- 
tions (solid curves). Nine ranges of halo mass from IO^'^Mq (top 
blue) to IO^^Mq (bottom red) are plotted. The dotted curves 
show the predictions given by integrating the mean M of our fit- 
ting formula (eq. 2). The lower four mass bins contain only haloes 
from the Millennium-II Simulation. For the upper five mass bins 
in which the haloes are drawn from both simulations, we use a 
solid circle to label the redshift above which only Millennium-II 
haloes contribute since the Millennium Simulation can no longer 
resolve haloes at such high z. The relatively smooth connection at 
the circle illustrates the consistency between the two simulations. 
Lower panel: Median mass assembly history M{z). For clarity, 
only two mass bins are plotted. Solid lines are from the simu- 
lations, dotted lines from the integration of the mean M from 
eq. (2), and dashed lines show the fits from Zhao et al. (2009). 

times greater at 2; = 2 (~ 2 x 10~* Mpc~^) than at a = 
(~ 3 X IQ-'^ Mpc-3). 



3.2.2 Mass Assembly Histories 

The top panel of Fig. 6 shows the mean mass assembly his- 
tory M{z) for nine bins of Mo (at z ^ 0) from 10^° to 
> 10^* (from top to bottom). The solid curves show the re- 
sults obtained from the main branch (i.e. the most massive 
progenitor) along the merger tree for all the z = haloes 
in the two Millennium simulations. The dotted curves show 
the M{z) obtained from integrating the fitting formula for 
the mean M in equation (2) from the present-day to some 



redshift z. The agreement is generally very good, in partic- 
ular at 2; < 8. A perfect agreement is not expected because 
the two quantities, M and M{z), are not determined from 
the same set of haloes in the simulations: the M statistics 
are obtained from all haloes of a given mass M at a given 
z, whereas the M{z) curves show only the mean mass of the 
most massive progenitors at redshift z for the z — haloes of 
mass M, which are a small subset of the haloes of the same 
mass that are present at z in the simulation boxes. In view of 
this difference, the agreement between the solid and dotted 
curves in the top panel of Fig. 6 is in fact quite remarkable. 
Over the large range of mass and redshift shown in Fig. 6, 
we have checked that the direct fits for the mean M{z) pro- 
posed in recent literature (e.g. Boylan-Kolchin et al. 2009; 
McBride et al. 2009) provide a good match at low z, but 
integrating the fit for {A4) in equation (2) provides a closer 
match at high 2, 

The solid curves in the bottom panel of Fig. 6 show the 
median, rather than the mean, mass assembly history ob- 
tained from the simulations for two mass bins centered at 
Mo = 10^'' Mq and IO^^Mq. We note that integrating our fit 
to the median M does not yield the median M{z) because 
unlike the mean M, the median and derivative operations do 
not commute. The median and mean M{z) are sufficiently 
similar, however, that we find integrating our mean M to 
yield relatively good agreement with the median M{z) (dot- 
ted curves). For comparison, the fit of Zhao et al. (2009) to 
the median M{z) is shown as dashed curves. Their fit ap- 
pears to be systematically lower than the Millennium results 
at 2 > 1. 

In principle, we can integrate the (mass-weighted) halo 
merger rate in equation (1) and obtain the portion of the 
dark matter accretion rate (M) in equation (2) that is due 
to mergers. As emphasized in Fakhouri & Ma (2010), how- 
ever, accretion of "diffuse" material (consisting of unresolved 
haloes and tidally stripped mass) also makes a non-negligible 
contribution to (M); equation (2) therefore can not be ob- 
tained solely from equation (1). 

3.3 Merger Statistics over a Halo's History 

In the last two sections we quantified the halo merger rates, 
the mass growth rates, and the assembly histories of haloes. 
These quantities can be combined to predict a number of 
additional useful merger statistics over a halo's history. 

3.3.1 Cumulative number of mergers 

One such statistic is Nm{(,min, Mo, zo, z), the total number 
of mergers that a halo of mass Mg at redshift zg has encoun- 
tered between 20 and an earlier 2. The mergers can be char- 
acterized by major or minor mergers by imposing a limit of 
^min on the progenitor mass ratio (evaluated at the redshift 
of the merger). These numbers are essential for making the- 
oretical predictions of galaxy properties that are impacted 
by mergers, e.g., the dynamics and stability of stellar disks, 
the star formation rate, and the color and morphology trans- 
formation due to mergers. 

Fig. 7 shows A^m(Cmin, Mq, zq , z) for the complete set of 
resolved haloes at 20 = (left), 1 (middle), and 2 (right) 
from the two Millennium simulations. In each panel, five 
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Figure 7. Mean number of mergers between redshifts zg and z experienced by a halo at zq = (left), 1 (middle), and 2 (right) from the 
joint dataset of the two Millennium simulations. In each panel, the solid and dashed curves represent mergers with a progenitor mass 
ratio (defined at the time of merger) of ^ ^ 0.3 and § ^ 0.1. For each mass ratio cutoff, five ranges of halo mass are shown (from bottom 
up): 10^° (blue), lO" (cyan), lO^^ (green), lO^^ (orange), and > (red). The lower -mass haloes are from the Millennium-II 

Simulation, whereas the cluster mass haloes are mainly from the Millennium Simulation. For mass bins in which the haloes are drawn 
from both simulations, we use a solid circle to label the redshift above which only Millennium-II haloes contribute since the Millennium 
can no longer resolve haloes at such high z. The fact that the curves connect quite smoothly are another indication of the consistency 
between the two simulations. 



ranges of Mq are plotted for redshift up to 12. Major mergers 
with ^ ^ 0.3 are shown by solid curves, while the more minor 
mergers with ^ J? 0.1 are shown in dashed curves. 

Fig. 7 shows that the mean trend of the number of merg- 
ers experienced over a halo's lifetime is a sensitive function 
of the halo mass and merger mass ratio, haloes of Milky- 
Way mass at the present day (green curves) have on aver- 
age experienced one major merger event (^ ^ 0.3) per halo 
since z ~ 2.3, and one merger with ^ ^ 0.1 per halo since 
z ~ 1. When extended to 2 ~ 7, these haloes have on aver- 
age encountered ~ 3 mergers with ^ ^ 0.3, and ~ 7 merg- 
ers with ^ J5 0.1. The formation redshifts as well as the last 
merger epoch for more massive haloes are both lower, a well- 
known trend in CDM-based cosmology (see, e.g., Lacey & 
Cole 1993, 1994). Cluster-sized haloes with Mq ~ IO^Mq, 
for instance, have on average experienced one major merger 
event since z ~ 1.2, and one merger with ^ ^ 0.1 since 
z ^ 0.6. 

It is possible to compute the cumulative number of 
mergers, Nmi^min, Mq, zo, z), shown in Fig. 7 from the fit- 
ting formula for the merger rate dNm/ d^/ dz in equation (1) 
and the mass accretion history M (2) obtained by integrat- 
ing (M) in equation (2). Specifically, these quantities are 
related by 

/■^ dN 
Nm{Un,Mo,Zo,z)^ dz di--^{M[z),i,z) . 

(3) 

Since we are interested in the number of mergers over a 
halo's past history, we must take into account the fact that 
a halo's mass generally decreases with increasing z, and that 
the merger rate depends on the halo mass (albeit weakly). 
The merger rate dN,n/d^/dz at redshift z in the integrand 
above therefore should be evaluated using the mean mass 
M{z) that a halo of mass Mo at 20 had at the earlier 2. 
The results are shown in Fig. 8, which is identical to the left 
panel of Fig. 7 except that we have added the theoretical 
curves (thin curves) for comparison. The agreement with 
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Figure 8. Same as the left panel of Fig. 7, with the addition of the 
predictions (thin curves) computed from eq. (3). The agreement 
with the simulation results (thick curves) is excellent, suggesting 
that eq. (3) can be used to make analytic predictions for merger 
statistics over a halo's history. 

the simulation results (thick curves) is excellent, suggesting 
that equation (3) can be used to make analytic predictions 
for merger statistics over a halo's history. 

3.3.2 Redshift of last major merger 

The redshift at which each curve in Fig. 7 crosses one merger 
event along the vertical axis is a useful quantity since it gives 
the mean redshift at which a halo has experienced its last 
major merger (LMM). The LMM redshift of a halo is closely 
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Figure 9. Cumulative distribution of the redsliift at which the last (i.e. most recent) major merger occurred in a halo's past history for 
haloes at zq = 0, 1, and 2 (left to right) in the two Millennium simulations. The vertical axis gives the probability that the last major 
merger occurred between zq and redshift z. The curves are defined the same way as in Fig. 7. 



related to its formation redshift and may be linked to the 
time at which the associated galaxy last experienced promi- 
nent star formation activity and morphological changes. To 
analyze this quantity further, we show in Fig. 9 the distri- 
bution of the LMM redshift for haloes at zq = 0, 1, and 2 
(from left to right). Within each panel, five halo masses and 
two ranges of ^ are plotted. The vertical axis gives the prob- 
ability that a halo at a given zq has had a last major merger 
between zq and z. 

Useful merger statistics can be read off from Fig. 9. 
For instance, 50% of present-day haloes have had a ma- 
jor merger (^ ^ 0.3) since z ~ 1, 1.8, and 3.4 for halo 
mass 10^'*, 10^^, and IO^^Mq, respectively. When more mi- 
nor mergers with ^ 0.1 are considered, the median redshift 
of the last merger is lowered to 0.4, 0.8, and 1.6 for the three 
masses. The assembly history of Milky Way-size haloes is of 
particular interest; see Boylan-Kolchin et al. (2009) and ref- 
erences therein for a detailed statistical study of this topic. 
For haloes of Mo ~ IO^^Mq today, the left panel of Fig. 9 
shows that ~ 31%, 53%, and 69% of these haloes have expe- 
rienced a major merger since z = 1, 2, and 3, respectively. 
For haloes of ~ lO^'^M© at 20 = 1 (middle panel), about 
50% of them have had a major merger since 2 ~ 2.7, and 
for haloes of the same mass at 20 = 2 (right panel), about 
50% of them have had a major merger since 2 « 3.7. 



4 SUMMARY AND CONCLUSIONS 

We have combined the halo catalogs from the two Millen- 
nium simulations to form an unprecedentedly large dataset 
for studying the merger statistics and assembly histories of 
dark matter haloes in the ACDM cosmology. The two sim- 
ulations provide, respectively, 11.3 x lO'^ haloes (between 
redshift and 6) and 7.5 x 10® haloes (between redshift 
and 15) above 1000 particles for our study. These haloes 
and their merger trees have allowed us to determine the dark 
matter halo merger rates and mass growth rates from 2 = 
to up to 2 = 15, for over five orders of magnitude in the de- 
scendant halo mass (10^° <Mo< W^^Mq) and progenitor 
mass ratio (10~^ C ^ !)■ For the small range of overlap- 



ping parameter space between the two simulations, we have 
found the agreement to be excellent. 

For the merger rates, the basic features reported in 
our earlier study based on the Millennium Simulation 
alone (Fakhouri & Ma 2008) are largely preserved in the 
Millennium-II Simulation. The mean merger rate per halo, 
dNm/d^/dz, is nearly independent of the descendant mass 
(Fig. 1 and 2) and scales as cx Mo°-^^^ at all redshifts. The 
merger rate in units of per redshift is nearly independent 
of redshift out to 2 ~ 15 (left panel of Fig. 3); the rate in 
units of per Gyr is therefore largely determined by the cos- 
mological factor of dt/dz and increases roughly as (1 -I- 2)^'^ 
at 2 > 1 (right panel of Fig. 3). Equation (1) provides an 
update on our simple analytical fitting form for the merger 
rate as a function of Mo, ^, and 2. 

For the mass growth rates of individual haloes, we have 
found the mean and median statistics (Fig. 4) to be well 
approximated by the simple fitting form of our earlier study 
(McBride et al. 2009). The updated coefficients based on 
the joint dataset from the two Millennium simulations are 
given by equation (2). The present-day mean and median 
rates at which a 10^^ Mq dark matter halo is accreting mass 
(at the virial radii) are 46.1 and 25.3 Mq yr~^, respectively. 
This rate increases nearly linearly with the halo mass, and 
increases with redshift approximately as (1 + z)^'"' at low 2 
and (1 -I- zf-^ at 2 > 1. 

We have also presented statistical quantities that track 
the merger histories of dark matter haloes cumulatively. 
Fig. 7 presents the number of major mergers experienced 
by haloes of various mass between redshift zg and 2 for 
20 = 0, 1, and 2. Fig. 9 presents the probability that a 
dark matter halo at redshift 20 will have last experienced 
a major merger at some earlier redshift 2. Much interesting 
and useful information regarding the contribution to halo 
growth made by major mergers can be read off these figures 
with ease. 

With the addition of results from the Millennium-II 
Simulation to our previous analysis of the Millennium Sim- 
ulation, the merger rate of dark matter haloes is now well- 
quantified for haloes with masses between 10^° and 10^^ Mq 
for redshifts 2 < 15, modulo the uncertainties inherent in 
halo definitions and in algorithms for handling fragmenta- 



© 2010 RAS, MNRAS 000, 1 13 



Merger Rates 11 



tion (see Appendix), for the cosmology used in the Millen- 
nium simulations. Several avenues remain open for future 
work, however. 

One obvious extension of the results in this paper is 
to consider the mergers of subhaloes themselves, as subhalo 
mergers can be more directly linked to galaxy mergers than 
can FOF halo mergers (Angulo et al. 2009; Wetzel et al. 
2009). Furthermore, the structure of the merger trees pro- 
duced for the Millennium simulations lends itself naturally 
to computing subhalo merger properties. While computing 
subhalo merger rates and connecting them to galaxy mergers 
presents additional challenges - in particular, the issues of 
assigning stellar masses to subhaloes, numerical resolution 
effects, and subhalo identification within larger FOF haloes 
- a thorough theoretical understanding of such rates is essen- 
tial for disentangling the relative contributions of merging 
and star formation to the growth of galaxies. 
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APPENDIX: COMPARISON OF DIFFERENT 

FOF MERGER TREES 

We refer the reader to Section 5 and Figure 8 of Fakhouri 
& Ma (2009) for a detailed discussion of the three basic 
operations - "snip," "stitch," and "snip" - that we have 
implemented and tested for handling the issue of halo frag- 
mentations during the construction of a merger tree for FOF 
haloes (see also Section 2.2 of this paper). Briefly, "snip" 
removes halo fragmentation events by severing the ances- 
tral link between the fragment subhalo and its progenitor, 
"stitch" places the fragment subhalo back into the FOF halo 
from which it emerged, whereas "split" removes the frag- 
ment subhalo's progenitor from its FOF halo, thereby gen- 
erating a new FOF at the progenitor redshift. 

Within the stitch and split algorithms, the operations 
can be applied either on a subset of fragments or on all frag- 
ments in a given FOF tree. Wc therefore subdivide each 
algorithm into two: stitch-oo vs stitch-3, and split-oo vs 
split-3. The stitch-oo and split-oo algorithms perform the 
given operation on all FOF fragments. This is done recur- 
sively from the redshift of fragmentation, going forward in 
redshift for stitch-oo, and backwards in redshift for split-oo, 
until there are no more fragments present in the simulation 
merger trees. As a result of this recursive process, stitch-oo 
identifies the first (highest-z) snapshot in which two sub- 
haloes join the same FOF to be their merger time, whereas 
split-oo selects the last (lowest-z) snapshot. 

These algorithms introduce some complications, how- 
ever. One particular problem faced by split-oo is the fact 
that there exists a firm cutoff' at z = 0, beyond which we 
do not have merger or fragmentation information. As a re- 
sult, although a fragment may actually finally merge beyond 
2 = 0, split-oo will incorrectly assign its final merger to an 
earlier redshift. This results in a pile-up of merger events 



at z = and, as we will show, artificially raises the low-z 
merger rate with respect to the high-z rate. Since there is 
no analogous hard limit at high z, stitch-oo does not suf- 
fer from this same behavior, and fragment mergers are re- 
distributed across all high redshifts evenly. On the other 
hand, any chance encounter between subhaloes that results 
in the FOF algorithm spuriously linking them together is in- 
terpreted as a real merger event by stitch-oo. The subhaloes, 
which may never interact again, arc nonetheless forced to 
join the same FOF group down to 2 = 0. 

The stitch-3 and split-3 algorithms are designed to limit 
the propagation effects of stitch-oo and split-oo. Stitch-3 per- 
forms the stitching operation on any FOF fragment that is 
observed to remerge with its progenitor FOF's main branch 
within 3 snapshots of the fragmentation event. Any frag- 
ments that do not satisfy this criterion are snipped, resulting 
in an orphan halo that may or may not later romorgc. The 
split-3 algorithm performs the split operation on any FOF 
fragment that is not a member of the main branch FOF 
at some point in the 3 snapshots before the fragmentation 
event. Again, fragments that do not satisfy this criterion are 
sruppcd. 

Neither stitch-3 nor split-3 adequately removes all re- 
merger events. Depending on the context this may be either 
a weakness or a strength: although the notion of halo re- 
mergers may be considered as multiple counting from a the- 
oretical perspective, observers will likely count as signatures 
all events that trigger mergers, regardless of whether they 
are the first or last entry. 

Moreover, both split-3 and stitch-3 have superior time 
convergence properties to the snip algorithm, in which the 
remergor problem is entirely unmitigated. Thus, stitch-3 and 
split-3 stand as intermediates between the snip and stitch- 
oo/split-oo algorithms. 

An immediate concern is whether the halo mass func- 
tion is heavily modified by the destruction/creation of FOFs 
due to the stitch/split operations. We have verified that 
these operations do not modify the mass function severely. 
For stitch-3 and split-3, the deviations are within 3% of the 
unprocessed mass function at all redshifts, while deviations 
of up to 10% exist for the stitch-oo (split-oo) algorithm at 
low (high) redshifts. 

Fig. 10 compares the five post-processing algorithms 
directly by presenting ratios of the per-halo merger rate, 
(INrn/d^l dz, as a function of progenitor mass ratio ^ at five 
redshifts (2 = 0.06, 0.5, 1, 2, and 4 from top bottom). Each 
column presents the ratio of the merger rate of a particu- 
lar post-processing algorithm (left to right: snip, stitch-3, 
stitch-oo, split-oo) to the merger rate extracted from the 
split-3 trees presented throughout this paper. Different mass 
bins are presented by different colored curves ranging from 
10^° Mg (blue) to 10^^ Mo (red). We note that though the 
region of overlap between Millennium (dashed) and Millen- 
nium 11 (solid) is small, there appears to be smooth continu- 
ation between these two sets of curves for all post-processing 
algorithms. 

The merger rates computed by all algorithms converge 
towards high z, though there is some residual disagreement 
with split-oo and snip at the ~ 20% level. There are, how- 
ever, distinct systematic differences among the algorithms 
when 2 < 4. Since we presented stitch-3 as our algorithm 
of choice for handling halo fragmentation in Fakhouri & Ma 
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Figure 10. Comparison of five algorithms used to handle halo fragmentation events in the FOF merger trees: snip, stitch-3, stitch-oo, 
split-oo, and split-3. Results presented throughout this paper are based on the split-3 tree. Plotted as a function of the progenitor mass 
ratio 5 (left to right) are the ratios of the merger rates, dNm/d^/dz, between each of the first four algorithms relative to the split-3 
results. Five redshifts are shown; z = 0,0.5, 1,2 and 4 (top to bottom). Within each panel, up to nine mass bins are shown: lO^'^Af0 
(blue) to > (red). The Millennium Simulation results are presented with dashed curves and Millennium II with solid curves. 

The systematic differences amongst the five algorithms are discussed in the text. 



(2008), we focus on the comparison of stitch-3 and split-3 in 
this section. The origins of the differences between the other 
algorithms and split-3 can be inferred from the discussion of 
the algorithms earlier in this section. 

The second column of Fig. 10 shows that stitch-3 and 
split-3 are in excellent agreement at all ^ for high mass haloes 
{Mo > 10^^). Low mass haloes, however, can show signifi- 
cant deviations in the merger rate. This is true especially in 



the major merger regime, where the merger rate predicted 
by stitch-3 is over 50% higher than split-3. This distinction 
was not detectable using Millennium alone, as the mass res- 
olution limited our analysis to Mo ^ IO^^Mq. 

To understand this deviation we have studied a subset 
of halo mergers in detail by analyzing halo tracks, velocities, 
and merger histories. In particular, we have constructed a 
number of criteria to determine whether a given merger is 
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actually a spurious encounter: if the relative velocity of the 
two haloes greatly exceeds the more massive halo's maxi- 
mum circular velocity, if the angle between the velocity vec- 
tors of the two haloes exceeds 70° at the time of merger, 
or if the FOF algorithm only associates the two haloes for 
two snapshots out of the eight snapshots centered on the 
merger snapshot, then the merger is deemed spurious. A 
qualitative look at three-dimensional halo trajectories finds 
that this criteria docs a good job of identifying chance halo 
encounters and premature mergers. 

For halo mergers with 1.1 x 10^° < Mq < 1.3 x 10^° 
and ^ > 0.1 at 2 = 0, we find that stitch-3 identifies 1,304 
mergers, while split-3 only identifies 738 mergers. Of these, 
505 mergers are in common, leaving split-3 with 233 merg- 
ers that are not in stitch-3 and stitch-3 with 799 mergers 
that are not in split-3. Of the 505 mergers in common, 
only 4 (0.8%) are deemed spurious by our criterion. Sim- 
ilarly, of the 233 mergers unique to split-3, only 12 (5.2%) 
are deemed spurious. Of the 799 mergers unique to stitch- 
3, however, 589 (73.7%) are deemed spurious. These spuri- 
ous mergers are primarily comprised of chance encounters 
in which the two otherwise unassociated haloes merge for 
a snapshot or two and then disconnect. While split-3 cor- 
rectly splits these events, stitch-3 does not and consequently 
inflates the merger rate. When these spurious mergers are 
removed, the remaining 210 mergers unique to the stitch- 
3 algorithm bring the stitch-3 and split-3 rates into close 
agreement. 

We note that depending on the context, one may choose 
one algorithm over another. Stitch-oo provides the first en- 
counter merger rate, but is known to link chance-encounter 
haloes that should not be linked. Split-oo provides the last 
encounter merger rate, but cannot be trusted for « < 1 and 
may incorrectly underpredict the merger rate due to spu- 
rious fragmentation. Split-3 stands in between both algo- 
rithms: it docs not propagate up and down the tree and 
does not heavily modify the distribution of FOFs, but it 
does double count some halo remerger events. This may be 
odious to the theorist, but may yield the most appropriate 
merger rate for comparison to observation. 
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